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The origin of predicted and observed anomalies in caloric 
curves of nuclei and other mesoscopic systems is investigated. 
It is shown that a straightforward thermodynamical treat- 
ment of an evaporating liquid drop leads to a backbending 
in the caloric curve and to negative specific heats in the two 
phase coexistence region. The cause is found not in the gen- 
eration of additional surface, but in the progressive reduction 
of the drop's radius, and surface, with evaporation. 



The thermodynamical equilibrium properties of first 
order phase transitions are completely describable in 
terms of the thermodynamic state variables associated 
with the individual separate phases. This is not the case 
in continuous phase transitions, where the two phases 
become progressively more similar as the critical point 
is approached. For this reason, in contrast with contin- 
uous phase transitions, first order phase transitions are 
"trivial," and interesting only in so far as they herald the 
appearance of a hitherto unknown or undescribed phase. 

Renewed attention to phase transitions has been gen- 
erated by studies of models with well defined Hamilto- 
nians with either short range interactions (e.g. the Ising 
model (|-[|] or the lattice gas model ]^-p^]) or incorporat- 
ing long range interactions such as gravitation or electro- 
magnetic interactions fuj.[u| Several of these stud- 
ies, microcanonical and canonical, were performed nu- 
merically, thus the results apply directly to finite (meso- 
scopic) systems. Features expected to disappear in the 
thermodynamic limit, if such a limit exists, were noticed 
and were claimed to be essential, characteristic indicators 
of phase transitions in mesoscopic systems |I^,|l3|,^0| ^9) . 
For instance, first order phase transitions were associated 
with anomalous convex intruders in the entropy versus 
energy curves, resulting in backbendings in the caloric 
curve, and in negative heat capacities p, p^ , ^3|,^0| - ^| . 

It is often claimed that these features appear only in 
microcanonical calculations and are thought to become 
lost or smeared out in the haze of canonical calculations 
^,^0-22 2^^(|. These anomalies have been attributed 
to a variety of causes, the foremost of which are surface 
effects (the energetic cost of an interface), and long range 
forces [g0[|4|-|2^]. Unfortunately however, the numerical 
nature of the calculations tends to make the identification 
of the causes of negative heat capacities rather problem- 
atic. 



In the experimental arena, a very elegant study has 
permitted the observation of the melting of clusters of 147 
sodium atoms and the associated negative heat capacity 

In the context of nuclear physics, microcanonical mod- 
els of nuclear multifragmentation have associated the 
anomalies of a convex intruder with the onset of multi- 
fragmentation pfj|^5| , p6| . Furthermore, lattice gas mod- 
els in the isobaric regime have also shown negative spe- 
cific heats in the coexistence region, where multifragmen- 
tation also appears (lj,[ll| ■ The question of whether the 
two transitions are related and possibly coincident with 
the liquid- vapor transition is still very much open. 

Recently, the claim has been made of an empirical ob- 
servation of these anomalies, such as negative heat capac- 
ities in nuclear systems |3l]] . These negative heat capac- 
ities have been inferred from the study of fluctuations in 
multifragmenting nuclear systems. Thus there is a great 
deal of interest in elucidating the origin of such anomalies 
in models as well as experiments. 

In particular, it would be highly desirable to ground 
any evidence for these anomalies, theoretical or other- 
wise, on thermodynamics itself, minimally modified to 
allow for the possible role of surface effects related to the 
finiteness of the system. 

In this paper we illustrate analytically how effects such 
as negative heat capacities can arise within a standard 
thermodynamic treatment. As an example, we consider 
the evaporation of a drop of ordinary liquid. Our only 
concern with mesoscopicity is the explicit treatment of 
the surface of the drop. 

Nuclear systems have long been associated with liq- 
uids, as testified by the success of the liquid drop model. 
In the spirit of the liquid drop model, the surface en- 
ergy introduces the simplest (and dominant!) correction 
to the bulk energy, leading to a model of one percent 
accuracy in systems as small as 40 nucleons and possi- 
bly smaller. A similar approach should hold for other 
kinds of clusters for which the surface energy may also 
be the first and dominant correction to the bulk prop- 
erties. Specifically, we will study the role of surface in 
generating the anomalies in the caloric curve. 

For many years, it has been known that a "mesoscopic" 
system, such as a tiny drop of liquid with radius r can 
be readily described in the pure thermodynamic limit 
p3|"p7| ■ The complete analogy between the liquid- vapor 
phase coexistence of a liquid in the bulk and for a drop 
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can be seen in Fig. [T], where the molar free energy F m 
at constant T is plotted versus molar volume V m for the 
two phases. The two free energy branches, for liquid 
and for vapor, can be considered completely indepen- 
dent. No interaction is assumed between the two phases. 
Coexistence becomes possible in the region of volume V 
where the overall free energy can be minimized through 
the common tangent construction. The equilibrium pres- 
sure of the saturated vapor is immediately given by 
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dF 
dV 



(1) 



No qualitative difference in the picture results by consid- 
ering a drop of finite radius r. The only difference is that 
the overall free energy of the drop (solid curves) as we 
shall see below, is higher than that of the bulk (dashed 
curves) and the equilibrium vapor pressure is correspond- 
ingly higher. 
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FIG. 1. A schematic plot of the molar free energy as a func- 
tion of molar volume for a liquid and vapor. Solid (dashed) 
lines demonstrate the behavior of a droplet (the bulk). Ver- 
tical lines show the molar volumes of the liquid and vapor at 
coexistence. Dotted curves represent superheated liquid and 
supersaturated vapor. 

The state of equilibrium between a liquid and its vapor 
can be described in the simplest way by the Clapeyron 
Equation 



dp _ AH m 
~dT ~ AV m T 



(2) 



where, p and T are the pressure and temperature, AH m 
is the molar enthalpy of vaporization and AV m is the 
difference of the molar volumes of vapor, V£, and liquid, 
yi 

Specialization to the case of a drop of radius r can be 
achieved by modifying the enthalpy to account for the 
surface energy pq] 
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where AH^ is the "bulk" molar enthalpy, S l m and Vj^ are 
the surface and volume of the drop and c s is the surface 
energy coefficient. 

Neglecting Vj^ compared to V£ and considering the 
vapor ideal, i.e. V£ = T/p, we can integrate Eq. (g), 
assuming also AH m to be constant. We obtain 



p = Po exp 



AH® 3c„Vl 



T 



rT 



P = Pbuik exp 



rT 



(4) 



(5) 



This equation contains all the thcrmodynamical infor- 
mation necessary to characterize the phase coexistence 
of the liquid drop of radius r with its vapor. The salient 
feature is the rise of the vapor pressure with decreas- 
ing radius. Fig. reffig:pt gives a map of the function 
p' = p'(T', r'), in terms of the scaled variables 
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For any given radius r, the function p — p{T, r) describes 
the equilibrium condition between the drop and its vapor. 
In other words, it is the phase diagram of the drop. The 
drop appears here as a "phase" defined by its radius r. 
A change in radius implies a change in phase. 




FIG. 2. The saturated vapor pressure as a function of the 
temperature for different droplet radii. The size of the open 
circles is proportional to the droplet radius. 

Let us now introduce some history and construct a 
caloric curve for a drop of radius tq. In order to operate 
at constant pressure po, we enclose the spherical drop 
in a deformable and expandable container to which we 
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apply an external pressure po. As the drop is heated, and 
before the vapor can appear, the temperature increases 
according to 



AH = ClAT, 



(J) 



where C l p is the liquid's heat capacity and is approxi- 
mately constant. When T reaches the value To at which 
the vapor pressure p(ro) = Po, the vapor first appears 
and it expands against the container. The heat of va- 
porization is absorbed at a rate H m (ro). However, as it 
evaporates, the drop sees its radius decreasing from its 
initial value, here chosen to be r' = 5.35. At constant 
temperature the vapor pressure would rise, but, at con- 
stant pressure, as we are now operating, the temperature 
decreases as shown in Fig. ||, as the system absorbs its 
heat of vaporization, so that, 
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After the drop has completely evaporated, the vapor can 
increase its temperature according to 



AH = ClAT 



(10) 



where C£ is the vapor heat capacity at constant pres- 
sure. The resulting caloric curve defined parametrically 
by Eq. (||) and Eq. Q and shown in Fig. |^ is rather 
interesting. It has a decreasing branch associated with 
the phase transition, along which the heat capacity is 
negative! (See Fig. ||.) 

As an aside we note that the scaled radius r' is just 
the ratio of the bulk energy to the surface energy. Thus 
for a nuclear system the range shown: 1 < r < 5.35, 
corresponds to a gold nucleus (A — 197) evaporating to 
a single nucleon. 

These rather extraordinary features are wholly due to 
the interesting but, in a way, accidental history of the 
decreasing radius with increasing evaporation. 

Because of the surface energy effect, each drop, of a 
given radius ro, is a separate phase in and of itself, dif- 
ferent from that associated with a different radius r\ . At 
fixed radius ro, nothing anomalous appears in Fig. |^; the 
pressure versus temperature curve, caloric curve and heat 
capacity are all perfectly normal. Anomalies arise when 
the system drifts from one radius to another, or from one 
phase to another. 

Typically, experiments p9|-|43| and calculations 
P, ^2p^j2tl| pT| ] heat a preassigned system with a cer- 
tain amount of energy AE or AH, and determine the 
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FIG. 3. The history dependent caloric curve of an evapo- 
rating drop at constant pressure. Dashed lines represent bulk 
behavior, solid line shows the drop's caloric curve. The scaled 
enthalpy is H' = H{AH° m /Zc a V^f /(An/Vi). 



resulting change in entropy S and temperature l/T = 
dS/dE\v', from the resulting caloric curve, phase coexis- 
tence diagrams are extracted. However, as shown above, 
the evolution of the system occurring during heating in- 
troduces complications in the construction of a phase di- 
agram from such a caloric curve. 

To avoid this problem an experimentalist or theo- 
rist would have to keep the radius constant or correct 
for its change while determining the vapor pressure as 
a function of the temperature, thereby eliminating the 
accidental aspects associated with the evolution of the 
system. The proper representation of all the thermo- 
physical properties associated with the coexistence liquid 
drop-vapor is that given in Eq. (^) and Fig. |^. 

The results obtained here are firmly grounded on ther- 
modynamics with a straightforward accounting of finitc- 
ness through the surface correction. They are exact in 
the limit in which the liquid drop model holds, namely, 
down to nuclei/clusters containing 20 or so constituents. 
They are completely general, as they do not depend on 
specific details of the system but rather on its gross prop- 
erties. In fact, they should be used as the paragon for 
lattice gas models and the like. In the limit in which these 
models represent liquid vapor coexistence, they must re- 
produce the present results. 

Even more importantly, this approach obviates the 
need for repeating numerical calculations for each indi- 
vidual system or drop size. All that is required is to 
determine the bulk energy (enthalpy) and the surface en- 
ergy coefficient of a give phase once and for all. 

It is not clear to us at the moment if the transition 
studied here and the anomalies associated with it have a 
direct counterpart in the allegedly observed phase tran- 
sition in nuclear systems |3l|]. It is however worth re- 
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FIG. 4. The history dependent heat capacity of an evapo- 
rating drop at constant pressure. Dashed lines represent bulk 
behavior, solid line shows the drop's caloric curve. The scaled 
heat capacity is C' p = C p (AH^/3c 3 V^) 3 /(iir/AH^Vi). 



peating that, once the constraint of constant pressure is 
enforced, the results described here are entirely general, 
as they apply to any small system undergoing solid- vapor 
or liquid vapor transitions. 

Anomalies in the heat capacities observed in mi- 
crocanonical calculations have been attributed to the 
increase in surface generated as additional liquid-vapor 
interface, e.g. in the formation of bubbles |?6). In the 
present case, however, and possibly generally, this con- 
clusion is not valid. In a finite system undergoing a 
liquid-vapor transition there is on average a decrease of 
surface as the evaporation proceeds. Any interior vapor 
bubble formation is disfavored by a Boltzmann factor 
exp(— AV/T), compared to the location of the same va- 
por on the outside of the drop, whose surface area ends 
up actually decreasing. Thus, the resulting anomalies 
are indeed surface related, but in a very different way. 

Several conclusions and caveats can be drawn from this 
treatment of an evaporating liquid drop: 

1. Mesoscopic systems can be dealt with within the 
context of standard thermodynamics, minimally 
modified to include the surface. 

2. Anomalous features such as a backbending caloric 
curve and attending negative heat capacities can be 
made to appear by allowing the system to "evolve" 
in parameter space (e.g. r). 

3. These features appear in a strict thermodynami- 
cal treatment and thus are not specific features of 
microcanonicity. 

4. These features are, in a way, accidental. They re- 
flect the evolution of the system in parameter space 



(here, the radius r). If the radius is kept fixed and 
the system is confined to a "single" liquid phase, 
the phase coexistence diagrams are completely or- 
dinary, and no new thermodynamics is evident. 
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